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A  HYBRID  MODEL  FOR  NEARSHORE  NONLINEAR  WAVE  EVOLUTION 

Kacey  L.  Edwards1  Jayaram  Veeramony2,  and  James  M.  Kaihatu3 

Field  and  laboratory  data  show  the  presence  of  an  equilibrium  in  the  high 
frequences  of  surf  zone  frequency  spectra  In  this  reglo„,  dissipation  of  £ 
requency  energy  equihbrates  with  nonlinear  interactions  that  transfer  energy 
pnmari  y  rom  low  frequencies  to  high  frequencies.  The  nonlinear  interactions 
that  describe  nearshore  wave  evolution  increase  model  complexity,  which  re- 
sulte  in  an  increase  in  computational  expense.  To  reduce  the  computational 
burden  of  a  nonlinear  wave  model,  it  is  combined  with  a  parameterization  for 
surf  zone  wave  spectra,  and  ifferent  implementations  of  the  parameterizations 
are  tested  and  the  results  show  that  the  predicted  spectra  and  Hrm.  compare 
»e  wit  ata  and  the  full  model  nonlinear  parameters  skewness  and  asym- 

sTggir  ^  aCCUrat6ly  PrediCted'  SeV8ral  We“  for  future  improvements  are 


INTRODUCTION 

It  is  well  known  that  as  waves  propagate  through  the  nearshore  en¬ 
vironment,  they  change  shape.  In  the  shoaling  zone,  simple  sinusoidal 

ZZTtT  mto  waves  wi,th  narrow’ peaked  crests  and  shallow’ broad 

f “V*™;  Furthermore,  near  the  area  of  breaking  the  waves  have  a  “saw- 

he  back  ofthe7’  ^  &0nt  ^  °fthe  Wave  is  and 

the  back  of  the  wave  js  more  gentle  sloping 

ener^Ir'Tr"”  "'arshor'  occur.  as  ,  result  of 

energy  transfers  between  different  frequency  components  of  the  wave  Al- 

s  from  7T  a11  frequencies- the  dominant  transfer 

move  r  Eh  ?TenCy’  In  addition-  dissipation  preferentially  re- 

transf  “7®’  fr°m  *7  hlgh  frequency  col*P°nents.  In  the  surf  zone,  the 
transfer  of  energy  to  high  frequencies  and  the  dissipation  of  energy  from 

higher  frequencies  reaches  an  equilibrium;  spectra  evolve  toward  a  broa“ 
featureless  shape  for  a  range  of  high  frequencies. 

The  transfer  of  energy  and  resulting  equilibrium  range  of  high  frequency 
aves  is  easily  seen  in  frequency  and  wavenumber  spectra.  As  depth  de- 
creases,  energy  at  the  peak  frequency  (/„)  or  peak  wavenumber  (kp)  of  a 

IprSSS  uST’  ^  R<SearCh  1009  Blvd..  Stennis  Space 

2 Jacobs  Engineering,  1009  Balch  Blvd.,  Stennis  Space  Center  MS  39529  USA 
Zachry  Department  of  Civil  Eiiffineerine.  A&K1  Tin  Ivnrcll-i  <  T  \  t  f  tt  ^  _  i 
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spectrum  decreases  and  the  tail  of  the  spectrum  obtains  a  gentler  slope. 
In  the  surf  zone,  the  slope  or  shape  of  the  spectral  tail  remains  relatively 
constant,  although  the  amount  of  energy  in  the  spectral  tail  may  decrease. 
The  transformation  of  spectral  shape  in  the  nearshore  environment  has 
been  observed  in  field  and  laboratory  data,  as  well  as  model  results  (Smith 
&  Vincent,  2003;  Berbers  &  Burton,  1997). 

Laboratory  Experiments 

Two  such  laboratory  experiments  are  the  experiments  of  Mase  k  Kirby 
(1992)  (referred  to  as  MK92)  and  Bowen  k  Kirby  (1994)  (referred  to  as 
BK94).  MK92  measured  random  waves  propagating  up  a  1  ^  20  bottom 
slope.  At  the  wavemaker,  the  Pierson-Moskowitz  spectrum  was  charac¬ 
terized  by  Hrm *  =  0,04  m  and  the  »  LG  Hz,  Therefore,  MK92  is 
an  example  of  mostly  deep  water  waves,  and  the  data  are  characterized 
by  weak  nonlinear  interactions  and  a  narrow  surf  zone.  BK94  measured 
random  waves  propagating  up  a  1  :  35  bottom  slope.  At  the  wavemaker, 
=  0*08  m  and  fp  =  0,225  Hz.  Therefore,  their  case  is  an  example  of 
mostly  shallow  water  waves.  Non! inear  interactions  have  more  of  an  effect 
on  the  transformation  of  wave  spectra,  and  the  surf  zone  is  wider,  relative 
to  MK92. 

Despite  weak  nonlinear  interactions  in  the  MK92  data,  changes  in  the 
spectral  shape  are  apparent  (Figure  1).  In  deeper  water,  where  nonlinear 
interactions  are  the  weakest,  there  is  little  change  in  the  spectral  shape. 
However,  as  depth  decreases  and  the  waves  shoal,  energy  at  the  low  fre¬ 
quencies  decreases,  and  energy  at  higher  frequencies  increases.  Therefore, 
the  slope  of  the  spectral  tail  decreases.  In  the  surf  zone,  significant  energy 
is  lost  from  the  spectral  peak,  and  the  spectrum  is  much  broader.  The  slope 
of  the  spectral  tail  is  flatter. 

Modeling  Nonlinear  Waves 

Because  waves  do  not  retain  a  simple,  sinusoidal  shape  throughout 
the  shoaling  and  surf  zones,  modeling  the  evolution  of  waves  through  the 
nearshore  environment  requires  a  more  complex  approach  than  linear  the¬ 
ory,  The  transfer  of  energy  that  causes  the  changes  in  the  shapes  of  waves 
and  spectra  is  facilitated  by  nonlinear,  triad  interactions.  These  three- 
wave  interactions  must  be  included  in  a  model  of  nearshore  wave  evolution 
to  accurately  predict  wave  shape  statistics,  which  in  turn  are  necessary  for 
sediment  transport  problems.  Unfortunately,  increasing  the  complexity  of  a 
nearshore  wave  model  by  including  triad  interactions  significantly  increases 
the  computational  expense  of  the  wave  model. 

To  reduce  the  computational  burden  of  a  nonlinear  wave  model,  we 
present  a  hybrid  model  that  combines  a  weakly  nonlinear,  dispersive  wave 
model  (Kaihatu  k  Kirby,  1995)  with  a  parameterization  of  the  surf  zone 
equilibrium  range  (Smith  k  Vincent,  2003),  In  addition  to  reducing  the 
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- - h  =  0.35  m 

- - h  =  0,15  m 

■  0  =  0,05  m 


1  2  3  4  5 

frequency  (Hz) 

S”"'“  ' 1  ~»«ct«d  during  ,hu  MK92  tub 

we  ““"P*"  «“>  *1“  hybrid  model’s  estimation  of 

sa.*r. 

HYBRID  MODEL 
Nonlinear  Wave  Model 

Many  types  of  models  are  capable  of  evolving  waves  through  the  nearshore 
environment,  and  each  model  has  its  adva„ta£s  and  diaadlta^ 
scnbed  herein  is  a  frequency  domain,  weakly  nonlinear,  dispersive  one- 

t—JT  m°.df  .(Kaihatu  &  Kirb*  1995)  (referred  to  as  NLMSE). 
Derivation  of  the  model  begins  with  the  boundary  value  problem  expressed 
in  terms  of  the  velocity  potential  (0)  in  only  one  horizontal  dimension  The 
governing  equation  and  boundary  conditions  are 

4>xx  +  4>zz  —  0; 

0!  =  —  hi0x; 

gn  +  4>t  +  I(0£  +  0*)  +  i,0It  =  0(e3); 

T)t  ~  4>Z  +  rjx0x  -  t 70tI  =  0(e3); 

where  0  is  the  velocity  potential,  e  (=  ka,  <£  1)  scales  the  nonlinear 


-h<z<  0, 

(i) 

*  =  -K 

(2) 

z  —  0, 

(3) 

£  =0, 

(4) 

the  nonlinear  i-prmc 
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tion,  and  the  subscripts  refer  to  derivatives  in  space  (x  and  z)  and  time  (f) 
respectively. 

The  boundary  value  problem,  governed  by  the  Laplace  equation  (1), 
describes  the  movement  of  an  inviscid,  irrotational  fluid.  Furthermore,  the 
problem  is  restricted  to  progressive  waves 


N 


n=l 


0  =  X  ~~fn^n(x,y)ei^lc"  dx_u'"t)  +C.C., 

—  *r 


where  N  is  the  number  of  frequency  components. 
The  evolution  equation  for  the  amplitudes  An  is 

jkccB)n 

'  o/j,  \  “n  t  Oxjdfj  — 


2  (kCCg)n 


rn-I 


mcc3)n 

N-n 


V  RAiAn-ie*  fik‘+^n-i-k„)dx 


+  2  X  SAtAn+iJfl** . 


[=1 


(5) 


(6) 


C  =  wnfkn  is  phase  speed,  and  Cg  =  d<jjnfdkn  is  group  velocity.  The  coef¬ 
ficients  ftandSarethe  subharmonic  and  super  harmonic  triad  interactions 
given  by 


*2 


+  {hi  -f  kn-t)(wn-ikt  +  tu*  Arn_f  )ujn] 


kn+l  +  (kn+i  -  ki )  (o?n+i  k(  +  uqfcn4*)u>n] 


y 


(7) 

(8) 


The  one  dimensional  model  is  initialized  offshore  with  a  wave  spec¬ 
trum.  It  moves  forward  in  space  determining  complex  amplitudes  for  a 
given  number  of  frequency  components  (N). 

Because  the  evolution  equation  is  derived  assuming  weak  nonlinearity 
and  mild  slope,  the  model  retains  the  linear  depth  dependence  and  linear 
dispersion  relation  given  by 


A- 


,  ,2 


cosh  kn(h  +  z) 
cosh  knh 


and 
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Aitnough  the  NLMSE  includes  terms  on  the  left  hand  side  that  combine 
shoahng,  refraction  and  diffraction,  it  does  not  account  for  energyTst To 

(SatInrto^egleI  href7eHK1hLtU  &  Kirb>'  (1"5)  add  the  third  term 
L  P  I  ,  0  ^  hand  Side  *he  evolution  equation.  The  dissipation 

term  is i  based  on  the  dissipation  model  of  Thornton  k  Guza  (1983)  They 
define  the  total  energy  dissipation  as  ^ 

Sy'rr  3/ 

*fc  =  ~WP9^jfHrms,  (11) 

where  B  and  7  are  free  parameters  with  values  1.0  and  0.6,  respectively 
and  fp  is  peak  frequency.  ’  pect,vely’ 

(lll'^SLiribuSiTf1  °'  ^  *  Kirby  ,1992)  th'  d“iP*«°«  model 
deni  component  “‘I”"-1'”*  depen- 


an“Qn"  +  (t) “*• 

(12) 

ano  =  Ft'b, 

(13) 

/P2X  i^i2 

Q«.  =  K  -  Gno]-*^ - , 

{14} 

Xa2  ki2 

n-1 

b  N 

(15) 

«Xc*w»ia 


,  ,  .  ,  —i-wuvj  vuui puiicuws,  ana  t  is  a  tree  parameter  that 

SffrZr  W  /erCe"tagA °f  thC  dissiPation  model  (11)  ^  distributed  to 
the  frequency  independent  function  (13). 

whirhhrlnrhtphapLd  Side°f  the  evolutio»  equation  contains  nonlinear  sums, 
hich  replicate  the  nonlinear  interactions  between  waves  of  different  fre- 

ES  m"d  r“te,y' ■“*  rnp,“ty  °f  ,he  £ 

mirnEfnT  7  eXpensive'  The  nonlinear  sums  require  a 

number  of  computations  on  the  order  of  N*.  In  addition,  the  NLMSE 

model  can  potentially  overpredict  energy  in  very  high  frequencies  (Figure 


.,w  .  ,h  °Ur  g°als  of  reducmg  the  computational  expense  of  modeling 

rZ  rlfrU  tlC  tlearShore  environment  and  improving  the  ability  to 
repjicate  high  frequency  energy,  we  employ  the  surf  zone  parameterization 
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Until  2003,  little  effort  was  put  into  parameterizing  surf  zone  waves. 
Following  the  development  of  deep  water  wave  parameterizations  (Resio, 
Pihl,  TYacy  k  Vincent,  2001)  and  previous  work  by  Thornton  (1977)  and 
Zakharov  (1999),  Smith  k  Vincent  (2003)  described  a  parameterization  for 
surf  zone(Jfrmj/fr  >  0.4)  wave  spectra. 

Relying  on  wavenumber  spectra  from  many  laboratory  and  field  ex¬ 
periments,  Smith  and  Vincent  (2003)  found  that  the  range  of  frequencies 
for  which  energy  transfer  and  dissipation  equilibrated  was  composed  of  two 
subranges.  They  call  the  two  subranges  the  Zakaharov  range  and  the  Toba 
range.  The  Zakharov  range  is  valid  for  low  wavenumbers  (kp  <  k  <  l/h)t 
and  the  Toba  range  is  valid  for  high  wavenumbers  {k  >  l/h).  Therefore! 
both  ranges  exist  only  if  kph «  L 

Each  subrange  has  its  own  shape  and  energy  level.  The  ranges  are 
defined  mathematically  as  power  laws  in  wavenumber  space  and  are  written 


h  =  0.35  m 


f  (a"1) 


Figure  2;  Results  from  the  NLMSE  model. 
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S(k)  = 

S(k)^0Zakk^- 

0T  ofca  =  OtTobah05 
0Zak  ~  CtZakk187 


kh  >  1 
kp  <  k  <  1/A 


(16) 

(17) 

(18) 
(19) 


a  s“br*nSes  intersect  at  k  =  1/A,  arpoba  =  oZak  =  0.0103  = 

•  addition,  the  dependence  on  depth  requires  that  the  energy  level 
decreases  as  depth  decreases.  gy 

Combining  the  Two  Components 

Use  of  the  parameterization  in  conjuction  with  the  NLMSE  model  in- 
res  a  reduction  in  computation  expense  because  the  number  of  frequency 
components  explicitly  determined  by  the  NLMSE  model  decreases  To 
combine  the  two,  use  the  relation 

$(fn)df  =  S(kn)dkn.  (20) 

r additi°n-  complex  amplitudes  are  determined  for  the  equilibrium 
range  using  the  relation 

A  -  l^(fn)df  i6 

4"-v~r~e  (2i) 

where  8n  is  a  random  phase  produced  by  a  random  number  generator. 
ALTERNATIVE  0roka 

Because  the  parameterization  was  developed  for  surf  zone  wave  spectra 
it  does  not  represent  offshore  wave  spectra  well.  At  offshore  locations,  the 
pe  of  the  parameterization  compares  reasonably  to  data;  however,  the  8 
coefficients  overestimate  the  energy  of  the  equilibrium  range  (Figure  3  ) 
Limiting  use  of  the  parameterization  to  the  surf  zone  region  works  well  but 

SaVingS-  Thk  is  especially  true  for  cases  similar 
to  MK92  where  the  surf  zone  is  relatively  narrow. 

An  alternative  approach  to  determine  the  beta  coefficient  was  explored 

™  was  ™tten  not  only  as  a  function  of  water  depth  but  also  as  a 

writing"  °f  thC  n°n  mearity  Parameter  Hrm,/h.  This  was  accomplished  by 

n.  —  pi  mj 

°  -  C~h~-  (22) 

C  is  a  constant  that  defines  a  as  a  fraction  of  the  nonlinearity  parameter. 

It  IS  given  by  the  relation 

0-ao7— —  (23) 

where  Qq,  ho,  and  0  represent  offshore  values,  and 
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BK94  data - $(K)Toba - S(k)z^ 


10 


0 


d  -  0.369  m 


k(m  ') 


Figure  3:  Comparison  of  Smith  &  Vincent  (2003)  parameterization  to  data  from  the  BK94 
experiment. 


&Toba,Q 


Mat 

/  S(k)Qdk 

Jk=:l/h _ 

h$*[N  k^2dk 

Jk^l/h 


(24) 


where  S(fc)0  is  the  wavenumber  spectrum  at  the  offshore  boundary* 

Frequency  spectra  resulting  from  the  NLMSE  model,  the  hybrid  model 
using  a  —  0*0103  (referred  to  as  SNV)  and  the  hybrid  model  using  (22) 
(referred  to  as  NALF)  are  compared  in  Figure  4.  As  expected,  offshore 
the  SNV  hybrid  model  provides  the  same  results  as  the  NLMSE  model.  In 
addition,  the  NALF  hybrid  model  compares  well;  the  results  overestimate 

^  £ - - 1  ^  _ _  XTT  % 


I  oun  t 
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hybrid  models  do.  However,  as  depth  decreases  the  NALF  hvhrid  a  i 

Ltss; Th:  snV  hybrid 

«.ergy.  Well  into  the  «,»,  both  hybrid 


h  =  0.369  m 


h  =  0.084  m 


4:  Frequency  „„„  renter*  ,ta  h(M  ^  —  ^  NLUSE  ^ 

SSJM:  TT °f  NLMSE  »Zut».,T2e 

T PT  W1'  ?  ^  d“*  “  *■  »“»«=  leeatioii 
tne  constant,  C,  depends  on  offshore  parameters  Thk  frsu*t\nn 

however,  .  not  constant  through  the  entire  donmin  CtI°n’ 

Comparisons  of  wave  parameters  determined  from  the  data  and  th* 

SSffSrf  TheectsNvath  \obn  r  the  frequency 
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tling  into  agreement  with  the  data.  The  NALF  hybrid  mode)  significantly 
underestimates  HTms  until  well  into  the  surf  zone. 


0  0.05  0.1  0.15  0.2  0.25  0.3  0.35  0.4  0.45 


Figure  5:  Wave  parameters  determined  from  NLMSE  and  hybrid  model  results.  Top:  //rTna, 
Middle:  skewness.  Bottom:  asymmetry. 


Comparisons  of  the  parameters  skewness  and  asymmetry  show  that 
neither  of  the  hybrid  models  reproduce  the  nonlinear  properties.  Accurate 
energy-  transfers  require  correct  phase  locking.  The  random  phases,  used 
to  determine  complex  amplitudes  in  the  equilibrium  range  do  not  ensure 
the  correct  phase  locking. 

PHASE  REPRESENTATION 

In  an  effort  to  improve  the  nonlinear  properties,  two  additional  phase 
representations  were  explored.  First,  the  offshore  phases  were  held  constant 
through  the  domain.  Then,  the  evolution  equation  was  used  to  evolve  the 
nhase  f.hrninrh  the  domain. 
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thatT°  °bUin  thC  eVOlUtiOD  equati°n  for  the  Phaa*s«  ***  began  by  assuming 

An  =ane'*"(*>.  (25) 

d,  S“b*titutine  <25)  into  (6).  we  solved  for  the  partial  derivative  of 

*n  and  used  the  result  to  specify  phase  for  the  equilibrium  range  ^thlt 


 1 

rn-1 

^1=1 

N-n 

+  2  ^  '  Sciian+ie*  $>,_ 

^  J 


(26) 


Figure  6  shows  skewness  and  asymmetry  results  from  the  SNV  hybrid 
e  using  random,  constant,  and  evolved  phases.  Neither  the  constant 

to  Soduce  theeV°1Vrd  PhaSe  aPP6arS  t0  impr°Ve  the  hybrid  model!s  ability 
to  produce  the  nonlinear  parameters  or  Hrma.  y 

SUMMARY 

.  JSSSLZTI  T”1”"™  thr°Ugh  and  s„f  ^  quires 

a  complex  model  that  is  computationally  expensive.  In  this  study,  the  com¬ 
plex  model  is  a  frequency  domain,  weakly  nonlinear,  dispersive  one  dimen 
smnal  wave  model  (Kaihatu  k  Kirby,  1995).  The  computational  expense  is 

“^lt7  thC  °f  frequenc-v  explicitly  deter- 

the  exctoded  fr°  lear  f  and  determinillg  the  complex  amplitudes  at 
the  excluded  frequencies  with  a  parameterization. 

The  surf  zone  parameterization  as  described  by  Smith  k  Vincent  (20031 

rhv3bLnlrithBthe  NTE  model  (Kaihatu  *  ***  mTtoZ 

3  model.  Because  the  parameterization  applies  to  only  the  surf 
one,  no  computational  savings  was  acheived.  However,  the  hybrid  model 
produced  a  better  replication  of  high  frequency  energy 

.  ,  lnf.an  3lteJrna*ive  d^cripton  of  0roba  allows  the  parameterization 

°U,5,d\,hC  rf““-  Therefore-  “>-P»tMk>„,l  burden 
evolving  waves  through  a  domain  is  lessened.  However  &r  t  's  denth 

(jendenre  d°ni1.inate*  the  Parameterization,  and  the  hybrid  model  does 
uot  compare  well  in  the  surf  zone* 

hvh  vi°del/fUltSTe  reiterated  in  the  computation  of  Hrm,.  The  SNV 
the  surfzom*  PH H'T  estimates  that  compare  well  offshore  and  inside 
Irl  f  H  r  °Werr’  here  15  a  region  of  overestimation  because  for  an 
a  of  the  domam,  the  SNV  hybrid  overestimates  the  low  frequency  energy 

For  most  nf  f-h*  *r*rr-  I  t  .  ,  .  4  J  energy. 
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Figure  6:  Wave  parameters  determined  from  the  SNV  hybrid  model  results  using  different 
phase  representations.  Top:  Hrm*.  Middle:  skewness.  Bottom:  asymmetry. 

Although  the  hybrid  models  give  reasonable  Hrma  values,  the  skewness 
and  asymmetry  are  clearly  inaccurate.  Using  a  random  phase  to  deter¬ 
mine  complex  amplitudes  in  the  equilibrium  range  interferes  with  the  phase 
locking  of  the  triad  interactions.  However,  neither  of  the  two  alternative 
approaches  of  determining  phase  for  the  parameterized  range  provided  any 
improvement  over  the  use  of  the  random  phase. 

Both  aspects  of  this  study  {the  energy  level  and  the  phase  represen¬ 
tation)  require  additional  investigation.  To  improve  the  0  coefficient,  the 
k  =  1/h  limit  will  be  explored  as  well  as  a  more  reliable  method  to  de¬ 
termine  the  constant  C  (23),  Furthermore,  perhaps  a  gentler  approach  to 
employing  the  parameterization  will  improve  the  results  and  reduce  the 
computational  expense.  Rather  than  forcing  the  spectrum  to  the  parame¬ 
terization,  maybe  it  would  be  better  to  guide  it  toward  the  parameteriztion. 
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